Description of the -0(3770) resonance interfering with the background 
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The parameters of the interfering i/)(3770) resonance should be determined from the data on the 
reactions e'*'e~ DD with the use of the models satisfying the elastic unitarity requirement. The 
selection of such models can be realized by comparing their predictions with the relevant data on 
the shape of the '(/'(3770) peak in the non-DD decay channels. Here, we illustrate this unitarity 
approach by the example of the most simple variant of the model of the mixed ?/>(3770) and ip(2S) 
resonances. When new high-statistics data become available, it will be interesting to test this clarity 
variant. 
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PACS numbers; 13.25.Gv, 13.40.Gp, 13.66.Jn 

In the recent paper [l| , we considered a few unitarized 
models available for phenomenological description of the 
e+e" — )■ DD reaction cross section in the V'lSTTO) reso- 
nance region. Such models allow us to avoid the spu- 
rious ambiguities in the interfering ?/;(3770) resonance 
parameters determination, which have been recently re- 
vealed by experimentalists when using unitarily uncor- 
rected parametrizations 

In this report we present the simplest working vari- 
ant of the model of the mixed -0(3770) and iJj{2S) reso- 
nances for the description of interference phenomena in 
the -i/'(3770) region. It was not discussed in Ref. 
Owing to own clarity and simplicity this variant can be 
tested, in the first place, in the treatment of new high- 
statistics data which can be expected from CLEO-c and 
BESIII H, [Mil on the V(3770) shape in e+e" ^ DD. 
Here we also concentrate great attention on the possi- 
bility of testing theoretical models by comparing their 
predictions with the relevant data on the shape of the 
-0(3770) peak in the non-DD decay channels, which are 
also expected from BESIII 

In constructing the model describing the process 
e+e^ — )■ DD, one must keep in mind that we investigate 
above all the D-meson isoscalar electromagnetic form fac- 
tor F]j. The phase of in the elastic region [i.e., be- 
tween the DD (w 3.739 GeV) and DD* 3.872 GeV) 
thresholds] is fixed by the unitarity condition equal to the 
phase Si of the strong P-wave DD scattering amplitude 
Ti in the channel with isospin / = 0, i.e.. 
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where J-'^^ and are the real functions of energy. A 
similar representation of the amplitude e+e^ DD used 
for the data description guarantees the unitarity require- 
ment on the model level jlj. The sum of the e+e~ — >■ DD 
reaction cross sections is given by 
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where s is the DD-pair invariant mass square, 

I'is) = [Po(s) +P+(s)]/^/i, Po,+(s) = ^ s/4 - m^o.+ and 

a — e^/47r = 1/137 (here we do not touch on the questions 
about the isospin symmetry breaking). Below, for short 
-0(3770) is denoted as ^p". 

Consider now the model which takes into account 
in F'^ and the contributions only from the ip" 
and -0(25") resonances. Owing to the common D^D^ 
and D^D^ coupled channels, the V" and the -0(25') 
can transform into each other (i.e., mix); for example, 
-0"— j'DD— >--0(2S'). The form factor F^, corresponding 
to the contribution of the mixed 0" and -0(25) reso- 
nances, can be represented in the following symmetric 
form P,[i3-[il: 
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where D^// (s) and D^(25) (s) are the inverse propagators 
of -0" and 0(25"), respectively. 



D^n(s) = mi,, - s - iy/sT 
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£^V(2S)(s) = ™^(2S) «\/sr^(2S)£)£)(s) > (5) 
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= 9ip{2S)-i[Dtl," {s)g^[2S)DD + n^"iA(2S)(s).9,/"'D£)] 

+9^"i[D^(2S){s)g^"DD + n^"^(2S')(s)5v-(2S)D5]- (8) 

The constants g^»jjQ, g^(^2S)DD, and g^»^, g^(2sh 
characterize couplings of the 0", 0(25) to the DD 
and virtual 7 quantum, respectively. The amplitude 
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n^"V'(2S)('S) describing the ■0" — '/'(2S') mixing has the 



form 
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Ren^"^(2S)(s) + « g;^j:^^ v{s). 

(9) 

Its imaginary part is due to the -0" DD ip{2S) tran- 
sitions via the real DD intermediate states. Substitut- 
ing Eqs. Q-dZl) and ^ into Eq. dH), it is easy to 
make certain that 'R-Dois) is a real function. Thus, the 
model can explain the dip observed in a(e~^e^ DD) 
near « 3.81 GeV (see Fig. [T]) by the zero in F^{s), 
caused by compensation between the ip" and ipiiS) con- 
tributions. Note that Ren^//^(2S)(s) cannot be strictly 
calculated. Its approximations, for example, by the ex- 
pression Co -|- sci, where cq and ci are free parameters, 
can be used as a resource for the fit improvement. Below, 
for simplicity we put Ren^//^(2S) (s) = 0. Then Eq. (|5]) 
takes the form 
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The curves in Fig. [T] correspond to 

9^jj"DD 

Eq. ®], g.^„^ 



3.794 GeV, 

±14.35 [i.e., r^„£,^(m|„) w 56.8 MeV, see 



- ip"DD\ 

= ±0.1234 GeV^ 



i.e., r 



'5^"7/(3"T.^") ~ 0.062 keV], and g^(2S)DD 
±20.11. In so doing, if 9^"'^9t[,"DD > (< 0), then 
9ij{2Sh94,{2S)DD < (> 0), see Eq. (HU]). The values 
'^i/'(2S) = 3.6861 GeV and 5^(23)7 = ±0.7262 GeV^ were 
fixed according the data [6|] and the relation r^(2s)e+e- = 

4^"'4(2S)7/(3"i^(2S)) = 2-35 keV. 

The values of the fitted parameters to^", 5^"dd; f-i^d 
g^".y can essentially depend on the model used for the 
description of the total contribution of the ijj" resonance 
and background. The analysis [1] indicates that the com- 
ponents of the e+e~ — >■ DD amplitude can be very differ- 
ent in the different models. For the model of the mixed 
-0" and ■0(25') resonances, the contributions of the com- 
ponents in question are shown in Fig. [1] by the dashed 
and dot-dashed curves. On the other hand, it is clear that 
the interference pattern in the ip" region depends on the 
reaction. Therefore, the selection of the theoretical mod- 
els should be carry out by comparing their predictions 
with the experimental data on the shape of the V'" peak 
for several different reactions. 

For example, after the fitting of the e+e^ -> DD data 
we all know about DD elastic scattering in the P-wave 
at the model level. 
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The corresponding cross section and phase are shown in 
Fig. [21 Unfortunately, these predictions are not possible 
to verify. However, there are other processes which can 
be measured experimentally. 
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Figure 1: The simplest variant of the model of the mixed tp" 
and ip{2S) resonances. The solid curve is the fit using Eqs. 
©-([Tni) to the data from BES 24, 25], CLEO 26], BABAR 
W^,iM, and Belle 29] for cr(e+e" -5> DO). The dashed and 
dot-dashed curves show the contributions to the cross section 
from the ■(/)" and tp{2S) production amplitudes proportional 
to the products of the coupling constants g^"-yg^ii on and 
5i/>(2S)75^(2S)DD, respectively; see Eqs. Q and Unj. For 
more details on the data see Ref. Q]. 




3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 
Vs (GeV) 

Figure 2: The predictions of the model with the 
mixed 1/)" and i/'(2<S') resonances. (a) The solid, 
dashed, and dot-dashed curves correspond to cr(Z)''_D° — 
D° D°)^3-k\ sin 5^{s)f/pl{s) and the -0" and ^(aS) contri- 
butions proportional to g^nj^f, and g'^(2S)DD ^1- CHi 
respectively, (b) The phase 5i(s). 
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We are interested in the interference phenomena in the 
tp" region in the reactions e+e~ — ^non-Z?l). We confine 
ourselves to the simplest non-DD final states, the form 
factors of which are determined by a single independent 
invariant amplitude. Such reactions are e"'"e~ — ^ 7XcOj 
777c, 777', J/^Tj, 077, and so on. 

The cross section for e'^e^^ab {ab = jXcO, 7^c, 7??', 
J /i^rj, (f)ri) in the t/j" region can be written as 



a{e'^e ah) — 



47ra^fcab(g) IP . ^|2 
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where kab{s) = ^/[s - {nia + rrib)'^] [s - {nia - mb)^] /{2^/s) 
and Fab{s) is the electromagnetic form factor of the ab 
system. Equation implies that the decay amplitude 
of the virtual timelike photon with the mass i/s into 
IXcO iXcO is the scalar meson) is given by 



(13) 



where £^{q) and el{k) are the polarization four- vectors 
of the intermediate (virtual) and final photons with four- 
momenta q {q^ = s) and fc, respectively; and, its decay 
amplitude into FO^ (0^ denotes a pseudoscalar meson 
and VQ~ — jrjc, jrj' , J/ipVi 4>'n) is given by 



aFvQ- (s) eM"^ffTe^(g)e^ {k)qakr . 
In the model under consideration we may write 



Fabis) 



D,pn{s)D^^^2S){s) - n^,"^(2S)(s) 
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where 



Tiabis) 

= 9iPi2Sh[F>i,"{s)gtl,(2S)ab + n^"^(25)(s)3-0"ab] 
+g4'"'y[F>4,{2S){s)gtlj"ab + n^"i/;(2S) (s)3i>(2S)ab] (16) 

and g^{2S)ab^ gip/'ab are the effective coupling constants 
of the ^{2S), ip" to the ab channel. These coupling con- 
stants are taken into account in Fab{s) in the first or- 
der of perturbation theory. Their relative smallness is 
caused by the electromagnetic interaction for the 7XcO 
and 777c channels, by the dynamics of the Okubo-Zweig- 
lizuka rule violation [15l - [l7| for the J/i/'f? and 4>r] chan- 
nels, and by a combination of the above reasons for the 
777' channel. 

As a first (rough) approximation, we suppose that the 
coupling constants for radiative transitions between char- 
monium states {cc)i — > 7 (cc) / [index i (/) labels ini- 
tial (final) state] and also those for hadronic transitions 
(cc)i — (cc) / h and radiative decays (ccji jh, probing 
the gluon content of light hadrons h, are real [Ol. IllL [l8l|. 
That is, we neglect the contributions of the real DD in- 
termediate states, taking into account which leads to the 
appearance of imaginary parts of effective coupling con- 
stants High-statistics studies of the e+e — >■ non- 
DD processes in the ip" region will show how this is jus- 
tified. Note that for the (cc)i — >■ 07; decay the DD loop 



rescattering mechanism (cc)i — > DD — )• is suppressed 
by the Okubo-Zweig-Iizuka rule. The phase of the (j)r] 
final state interaction is unknown. However, this phase 
is common for different contributions to e^e~ — >■ (j)r] and 
does not appear in the cross section. At this stage, we 
do not take into account the interference between the 
e+e^ — >■ (cc) — >■ 077 amplitude and the background from 
the light quark production e+e~ — )• (ss) — > ^77. With 
the above assumptions, the effective coupling constants 
9ip{2S)4>7] and g^"(i,ri will be considered to be real as well. 

Table I presents information about the '0(25') Q and 
ip" [1, [19I - I22I resonances in the ab decay channels, which 
we use to construct the corresponding mass spectra. The 
values for g^(2S)ab indicated in the table are obtained, 
up to the sign, from the data on the ip{2S) ^ ab decay 
widths by the formula 



_ 9il,(2S)ab ,3/2 N 
V(2S)a6 - l^ab\'^^(2S)) 



(17) 



which implies that the amplitudes of the '0(25') — )■ 7XcO 
and ^(25) VQ^ decays have the form 

5^^(25)7X00 4^^'^'^\,q-hZ{k){q ■ kg,,,, - k^q„) 

and 5-0(25) vo- £f.,^aTet^^^'^\q)e^ {k)q^kr , 

respectively. The relative signs of the constants g^{2S)ab 
and g^"ab are unknown. Therefore, the relative signs 
between the first and subsequent three terms in Eq. (|16p 
(they are controlled by signs of the coupling constant 

products) can be chosen in two ways: (H h) or ( — h 

— ). Here, we took into account the above-mentioned sign 
correlation between g^"^g^„DD and 5v>(2S)75v(2S)dd- 

The existing information about the 0"— >7XcOi IVci 
777', J/'ipr], (j)rj decays are very poor. The CLEO Col- 
laboration measured the reactions e+e^— ?>7XcO (19i] . 
e+e~ — ?► J/V"?? dH: and e+e~ 4>ri [l^l at a single point 
in energy ^/s = 3773 MeV (at the supposed maxi- 
mum of cross sections). The approximate values for 
(T(e+e~ — )■ ab) are presented in Table I and Fig. [3] by the 
points with the error bars. They allow us to roughly es- 
timate the coupling constants 5v"7Xco ~ ±0.608 GeV~^, 
g^nj/^^ « ±0.0375 GeV-i, g^n^^ « ±1.1 x IQ-^ Gc\-^ 

for case (-1 h) aiid g^n^^^^ sa ±0.721 GeV^"'^, g^"j/^rj ~ 

±0.065 GeV^i, 5^"^^ « ±1.11 x 10"^ GeV"! for case 
(- + -), by using Eqs. ([H]), 1^^, and and con- 

struct the corresponding cross sections as functions of 
energy. 

The solid curves in Figs. [3] show the cross sections for 



■IXcO, 



■ J/tprj, and e 



■ (j)T]; the dashed 



and dotted curves show the contributions from the 0" 
and 0(25) resonances proportional to [see Eq. (HI 



[5i/;"7£'0(2S)(s) + 50(2S)7n^"V(2S)(s)]5^A"Qb 

and [5v(2S)7-D0" (s) + 5^"7H^//^(2S) (s)]5v'(2S)ab , 

respectively. The values of each of these contributions 
to Fab{s) change from reaction to reaction according to 
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Table I: Information about the il){2S) fil] and ip" [T^ - [2^ resonances in non-DD decay channels (ab). 



ab 


B{^{2S) ah) 


r^(2S)ab (keV) 


g^(2S)ab (GeV ^) 


B{iP" -> ab) 


r^„ab (keV) 


a{e^e ah) (pb) 


7XcO 


(9.68 ±0.31)% 


29.4 ±0.9 


±(0.250 ± 0.004) 


(7.3 ±0.9) X 10-3 


172 ± 30 


72 ±9 




(3.4 ±0.5) X 10-=^ 


1.03 ±0.15 


±(1.22 ±0.09) X 10-^ 








W 


(1.23 ±0.06) X 10-^ 


0.374 ±0.018 


±(1.67 ±0.04) X 10-3 


< 1.8 X 10-'' 






J-tprj 


(3.28 ±0.07)% 


9.97 ±0.21 


±(0.218 ±0.002) 


(9 ± 4) X 10-"* 


21 ± 10 


5.5 ±2.5 


<Pv 


(2.8tJ:«) X 10-5 


(8.5l^'0) X 10-3 


±(2.7l»;5) X 10-* 


(3.1 ±0.7) X lO-* 


7.4 ± 1.6 


2.4 ±0.6 





3.74 3.76 3.78 3.8 3.82 3.84 3.86 3.£ 

Vi" (GeV) 



3.74 3.76 3.78 3.8 3.82 3.84 3.£ 
\^ (GeV) 




3.74 3.76 3.78 3.8 3.82 3.84 3. 
Vi^ (GeV) 



3.74 3.76 3.78 3.8 3.82 3.84 3.8 
Vi" (GeV) 



S 3 




3.74 3.76 3.78 3.8 3.82 3.84 3.86 3.88 
4% (GeV) 




3.74 3.76 3.78 3.8 3.82 3.84 3.86 3.£ 
•fi (GeV) 



Figure 3: The cross sections for e^e ^7XcO, e.'^e -^J/ijiri, 
and e'^e~^(j>rj (left) for case (H h) and (right) for case 

(-±-)- 



changes of g^"ab and g^(2S)ab- At the same time, their 
s-dependence does not change, as it has aheady been 
determined by the model parameters found from fitting 
the e^e~ ^ DD cross section (simultaneous fits to the 
data on the reactions e+e" — > DD and e+e~ — >■ non-Di? 



is yet to come). Note that the cross section for e+e- (t>rj 
is completely dominated by the V" contribution. Note 
also that the Belle Collaboration has recently measured 
the cross section for e+e" — ?> J H'ri between -^s — 3.8 GeV 
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Figure 4; The cross sections for e'^e -^frjc (left) and 
e+e"— >777' (right). 



and 5.3 GeV [23|. Unfortunately, the data for 3.8 GeV< 
y/s < 4 GeV have large errors, which does not allow us 
to extract any useful information. 

The cross sections for e+e- — >7?7c and e+g- -^^rj' in 
the ip" region are unknown. Using information about 
the ?A(2S') from Table I, we estimate the cross sections 
at g^'i-yrjc =9ip"'yri' =0. The results are shown in Fig. 0] 
by the dotted curves. Here, as in the case of the dotted 
curves in Fig. [3J the resonant enhancement on the tails 
of the 'ip{2S) contribution arises owing to the -ip" — V'(2'S') 
mixing. If we put T^/i^j^^ « 1 keV [l^, which corresponds 
to g^"-yr)a ~ ±1 X 10~^ GeV~^, then (T(e+e~ — >■ 777c) takes 
the form shown in the left plot in Fig. |4] by the solid 
curves for cases (H h) and ( — \ — ). 

The above examples tell us that the mass spectra in 
the '0" region in the non-DD channels can be very di- 
verse. Therefore, we should expect that the data on such 
spectra, together with the e+e~ — >■ DD data, will impose 
severe restrictions on the constructed dynamical models 
for the ip" resonance interfering with the background. 
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